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Abstract 

WL study the dynamics of a plasma of charged relativistic fermions at very high tem¬ 
perature T ^ m, where m is the fermion mass, coupled to the electromagnetic field. 
In particular, we derive a magneto-hydrodynamical description of the evolution of such a 
plasma. Bh show that, as compared to conventional MHDfor a plasma of non-relativistic 
particles, the hydrodynamical description of the relativistic plasma involves new degrees 
of freedom described by a pseudo-scalar field originating in a local asymmetry in the den¬ 
sities of left-handed and right-handed fermions. This field can be interpreted as an effec¬ 
tive axion field. Taking into account the chiral anomaly we present dynamical equations 
for the evolution of this field, as well as of other fields appearing in the MHD description 
of the plasma. Due to its non-linear coupling to helical magnetic fields, the axion field 
significantly affects the dynamics of a magnetized plasma and can give rise to a novel 
type of inverse cascade. 


1 Introduction 

It has been suggested [1,2], more than a decade ago, that the chiral anomaly may be accom¬ 
panied by the appearance of some new degrees of freedom described by an effective axion 
field that couples to the index density of the electromagnetic field. More recently, it has been 
understood [3] that this axion field plays an important role in the dynamics of the electromag¬ 
netic field in a relativistic plasma of charged fermions, (e.g., such as a plasma described by 
the Standard Model at high temperature and/or large matter density). In the spatially homoge¬ 
neous case, the time derivative of the axion field is proportional to the difference between the 
chemical potentials (densities) of left-handed and right-handed charged particles [1,2]. As a 
consequence of the chiral anomaly, a left-right asymmetry, as described by a non-vanishing 
time-dependent axion field, affects the dynamics of the electromagnetic field (see also [4-7]) 
and gives rise to an instability of solutions of the coupled Maxwell-axion equations [1]. 

This kind of instability is encountered in different physical systems exhibiting a left-right 
asymmetry. For example, it has recently also appeared in the electrodynamics of topological 
insulators; (see, e.g., [8]). Instabilities leading to the growth of helical magnetic fields in the 
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plasma of relativistic charged particles, as discussed in this paper, are expected to have played 
an important role in the dynamics of the early Universe, assuming that a potent source of chiral 
asymmetry was present; (see, e.g., [1, 3, 5, 9-1 1]). Magnetic fields pervading the primordial 
plasma at early times would have had significant effects on many processes, such as baryo- 
genesis [12], primordial nucleosynthesis [13, 14], CMB physics (see [15] for a review), the 
growth of various inhomogeneities (see [16], and refs, therein), etc. Thus, understanding the 
dynamics of primordial magnetic fields in the presence of sources of left-right asymmetries, 
as described by time-dependent effective axion fields, is essential in attempts to describe the 
evolution of the early Universe. 

An important question in cosmology is whether traces of primordial magnetic fields could 
have served as seeds for galactic magnetic fields [17, 18] or given rise to the large-scale mag¬ 
netic fields observed in voids [19-21]. The astrophysical origin of the magnetic fields in 
voids is still under active investigation; (see e.g. [22, 23]). If their origin turned out to be 
primordial then the difficulty that causally independent homogeneous magnetic fields survive 
[24-26] would have to be faced. To cope with this difficulty one might argue that the mecha¬ 
nism of generation of those fields is “super-horizon”, presumably related to inflation [27, 28]. 
However, the generation and survival of somewhat sizeable magnetic fields during inflation is 
difficult to reconcile with slow-roll conditions [29-32]. 

In an expanding Universe filled with Standard-Model particles, chirality flipping reac¬ 
tions, which erase chiral asymmetries, are in thermal equilibrium at temperatures below ~ 80 
TeV [33]. It has therefore been expected that an initial growth of magnetic fields, as derived 
from the chiral anomaly, could only have taken place at such truly enormous temperatures (if 
they ever existed in the early Universe). In an expanding Universe, the wavelength of any 
perturbation grows like the scale factor a{t) oc (with a < 1, for any non-inflationary expan¬ 
sion), and the distance, Rh, to the horizon in the Universe increases linearly in time, Rh oc t. 
As a result, even if a magnetic field were generated during some very early epoch, with corre¬ 
lations extending over a horizon-size scale (maximal scale possible for any causal generation 
mechanism), it would soon become stronger at sub-horizon scales. But it is known that sub- 
horizon-scale magnetic fields die out because of dissipation caused by resistivity and viscosity 
effects [24, 34-38]. It has been argued that the only way for causally generated magnetic 
fields to survive is the turbulence-driven inverse cascade [24, 25, 39, 40] of helical magnetic 
fields - a process that makes the characteristic scale of the magnetic fields grow faster than 
red-shifting does. At present, the efficiency of this type of mechanism is actively discussed. 

We have demonstrated previously [3] that there exists another (new) mechanism for the 
transfer of magnetic helicity from short wavelengths to long wavelengths. This mechanism is 
not related to magneto-hydrodynamical turbulence, but has its root in the chiral anomaly. We 
have then argued that, because of their coupling to the effective axion field (which describes 
chiral asymmetries), helical magnetic fields may survive down to rather low temperatures of 
a few MeV, even though chirality flipping reactions (caused by mass terms for the charged 
particles) are in thermal equilibrium. It has also been argued in [41] that, as a consequence of 
the parity-violating nature of weak interactions and the chiral anomaly, quantum corrections in 
states of finite lepton- or baryon number density may give rise to an asymmetry between left- 
handed and right-handed particles. A time-dependent effective axion field is then generated, 
which, in turn, excites a magnetic field [1]. One concludes that a homogeneous, isotropic 
stationary state of the Standard-Model plasma might be unstable. Thus, one has to study 
dynamical features of this system, in order to identify its true equilibrium state. 
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Large-scale electromagnetic fields in the plasma excite macroscopic flows of matter con¬ 
sisting of charged particles. The dynamics of low-energy modes of the plasma can be de¬ 
scribed with the help of the equations of relativistic magneto-hydrodynamics (RMHD). One 
might expect that the hydrodynamical description of a magnetized relativistic plasma does not 
differ in an essential way from the one of its non-relativistic cousin, the only important differ¬ 
ence being a different relation between energy density, p, and pressure, p; ( for reviews see, 
e.g., [24, 42, 43]). However, the chiral anomaly actually leads to an important modification 
of the Navier-Stokes equation; see e.g. [44-46]. Furthermore, helical electromagnetic fields 
affect the evolution of the difference of chemical potentials of left- and right-handed particles, 
which becomes (space- and) time-dependent [3]. As a result, not only the Navier-Stokes equa¬ 
tions, but also the connection between the electric current density and the magnetic field are 
modified in accordance with the chiral anomaly, and an equation of motion for an effective 
axion field must be added. For a spatially homogeneous state, the time derivative of this axion 
field is nothing but the chiral chemical potential, (i.e., the difference of chemical potentials for 
left- and right-handed particles). 

A magneto-hydrodynamical description of a chiral relativistic plasma of charged massive 
particles, as studied in this paper, provides an approximate description of the Standard-Model 
plasma at high temperatures and positive matter density; (additional effects not related to the 
chiral anomaly must, however, be taken into account, too). We expect that our magneto- 
hydrodynamical description of the plasma is a useful and reliable tool to determine properties 
of (local) thermal equilibrium of the plasma at non-vanishing lepton- and baryon densities. 

Apart from applications in studies of the Early Universe, our magneto-hydrodynamical 
approach may lead to a reasonably accurate description of astrophysical systems, such as rela¬ 
tivistic jets, and of heavy ion collisions [45] . Moreover, it may change our current views of the 
origin and evolution of primordial magnetic fields and their relationship to the presently ob¬ 
served large-scale cosmic magnetic fields [19-21]. Our results and the effect reported in [41] 
could affect our current understanding of cold but very dense systems, such as neutron stars. 
Last but not least, our magneto-hydrodynamical approach may be useful in the analysis of 
models of beyond-the-Standard-Model particle physics, such as the z/MSM [47-49]. 

We stress that the effects studied in this paper do not involve any “new physics”; they 
can be described within the Standard Model. The axion field introduced below might not 
be a fundamental field but, rather, an emergent one that appears in the description of certain 
states of the Standard-Model plasma. If instead a fundamental axion does exist, then the high 
temperature states of Standard Model plasma augmented by this new field (as in [1, 2]) can 
also be described by the chiral RMHD and the effects, discussed in this paper could emerge in 
a wider class of states. 

1.1 Main goals of the paper 

The physical system of primary interest in this paper is a relativistic plasma of very light 
charged particles at a temperature T large enough that the masses of the lightest charged par¬ 
ticles, the electrons and positrons, can be neglected, (i.e., set to 0). For simplicity, one may 
consider a plasma consisting only of electrons and positrons, with a charge-neutralizing back¬ 
ground of protons (and neutrons). But our analysis readily extends to more general systems. 

Let J^{x) denote the usual axial vector current density, where x = (f,x) are coordinates 
of a space-time point; E{x) and B{x) denote the electric field and magnetic induction, respec- 
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tively. (In our notations, we will not distinguish the quantized eleetromagnetie field from its 
expeetation value inphysieally interesting states.) Furthemore, A(x) = (Ao(x) = (j){x),A{x)) 
is the electromagnetic gauge field, with 0 the electrostatic potential and A the vector potential. 
The chiral anomaly says that is not a conserved current; more precisely that 

2(y ^ ^ 

df^J^ix) = —E ■ B{x) , (1) 

TT 

where a is the fine structure constant. Introducing the current density 

^ ( 2 ) 

where Fp^ are the components of the electromagnetic field tensor, we find that 

= 0, (3) 

i.e. that is a conserved current density. Unfortunately, it is not invariant under electromag¬ 
netic gauge transformations, A^ ^ A^ — dpT], where r] is an arbitrary function on space-time. 
However, the charge 

Qs = J d^xJ^{x,t) (4) 

t=const 

is not only conserved but also gauge-invariant. In order to describe thermal equilibrium states 
of the plasma, we are advised to introduce an ‘"axial chemical potential”, conjugate to Q^. 
The chemical potential conjugate to the conserved total electric charge is tuned in such a way 
that the system is neutral, (i.e., the charge density of electrons and positrons in the system is 
canceled, in average, by the one of the heavy charged hadrons, in particular the protons). 

We introduce the pseudo-scalar density 

P5{x,t) = {J°{x,t))T,^, 


and its spatial average 

gs = p5{x,t), 


where we denote by {{■))t,ii 5 the equilibrium state of the system at temperature T and axial 
chemical potential p^, and by (■) we indicate spatial averaging. As one would expect, there is a 
response equation relating the average axial charge density, q^, to the axial chemical potential 
/X 5 of the form 


q^{T,p^) ^ /is 


dq^ 

dp5 


M5=0 


6 


(5) 


with q^{T, /is = 0) = 0. A “derivation” of this relation can be found in cosmology text books.* 
Using the chiral anomaly, we see that Eq. (5) implies the equation 


/is = |^?5(T,/i5) = 


( 6 ) 


a relation derived and explained in [3]. (The dot indicates a derivative w.r. to time, and we 
neglect terms ~ T!) It has been shown in [1, 7] that, in the presence of a non-vanishing 

'it is assumed that in processes such as e’*' -|- e“ V-)- 7 -I- 7 and -|- 7 <->■ -I- 7 the rates of direct and 

inverse reactions are identical, see, e.g., [50], Sect. 3.4 and Chapter 5. 
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magnetic induction B{x), a non-zero axial chemical potential, ^ 0, induces an electric 
current density given by so that the complete expression for the electric current 

density, j, of the plasma is given by 


a 


j{x) = (J(a;))T ,/.5 = -l^bB{x) + aE{x) , 

TT 


(7) 


where a is the Ohmic conductivity.^ In [1, 7], the term ^fi^B on the right side of (7) has been 
derived from the chiral anomaly using current algebra; (see [51]). 

Inserting (7) in Ampere’s law and applying Faradey’s induction law, one readily finds the 
equation^ 


ot a 71 a 


( 8 ) 


As shown in [ 1 , 5], the term proportional to /xs in Eq. ( 8 ) is responsible for an instability of the 
solutions of ( 8 ), namely an exponential growth of the magnetic induction B, which leads to a 
non-vanishing density E-B, (i.e., to a non-vanishing helicity of the electromagnetic field). By 
eq. ( 6 ), this, in turn, leads to a temporal variation of the axial chemical potential /is and, hence, 
to a change in time of the pseudo-scalar density gs that describes the asymmetry between the 
charge densities of left-handed and right-handed charged particles. (We note that the electric 
field is damped by Ohmic losses, which affects the time dependence of E ■ B.) 

Equation ( 6 ) implies that, in general, /is depends on time. The anomaly equation (1) shows 
that the chiral density P5(T', t) usually not only depends on time, but also on space, x. It is 
therefore necessary to generalize relation (5) to one where /is not only depends on t, but also 
on space x; to then derive the correct form of the Maxwell equations in the presence of a space- 
and time-dependent axial chemical potential, /is(a:, t), and to derive the equations of motion 
for this field. To find such a generalization is one of the main goals of this paper. We will 
also study how to couple our system of equations to the equations of motion of a relativistic 
fluid, thereby deriving the correct equations of relativistic magneto-hydrodynamics (RMHD). 
In Sect. 4 we discuss some simple solutions of these equations and sketch possible applications 
to plasma physics, astrophysics and cosmology. 

As was pointed out already in [1], electrodynamics coupled to an axion field, 9^, can 
be used to describe the dynamics of a homogeneous system (with 6*5 = 6*5 (f) independent 
of x) in the presence of an axial chemical potential p^{f) = 9^{t) satisfying Eq. ( 6 ). In 
this work, we show that, for an inhomogeneous chemical potential p^{x,t) (defined below. 
Sect. 2), equations analogous to Eqs. ( 6 ) and ( 8 ) can be written as local field equations for the 
electromagnetic field and for an effective axion field 9^{x), with 9^{x,f) = p^{x,t), (rather 
than in terms of the axial chemical potential p^{x,f)), as already pointed out in [1]. The 
axion field 6^{x, t) interacts with the electromagnetic field in the form of a term O^E ■ B in 
the Eagrangian density, but its dynamics is not necessarily described by a relativistic wave 
equation, but might be governed by an inhomogeneous diffusion equation. 

So far, the effective axion field 6*5 has been treated as a classical field. The role played 
be the axial chemical potential p^{x,f) in the description of inhomogeneous local thermal 
equilibrium states suggests that and hence the axion field 6^ does not describe hamiltonian 
degrees of freedom that have to be quantized; but that they should be understood as c-number 


^For simplicity, the state of the plasma is assumed to be isotropic and homogeneous. 

^The standard approximation made in MHD is to neglect Maxwell’s displacement current and hence the 
second derivative, on the right side of (8); see, e.g., [52]. 
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fields labeling a family of generally inhomogeneous non-stationary states of the quantum- 
meehanieal plasma - as it is the ease of the loeal temperature (field) T{x). 

One eould also eonsider a system with an axion field does deseribing dynamieal (hamilto- 
nian) degrees of freedom that must be quantized. In order to shed new light on this important 
issue, we show, in Appendix A, how an effective axion field emerges in theories with extra 
dimensions. In particular, we consider a model of (quantum) electrodynamics on a slab in 
five-dimensional Minkowski space in the presence of a 5D Chern-Simons term. In this theory, 
the axion appears as a component of the electromagnetic gauge field (in the fifth direction 
transversal to the boundaries of the slab), and the state function T(a;) has an interpretation as 
the inverse width of the slab and is therefore related to the geometry of five-dimensional space- 
time; (see also [53]). Comparing this 5D formulation with the analysis presented in Sects. 2 
and 3, one may wonder whether, in the end, space-time geometry might merely encode prop¬ 
erties of the state of quantum-mechanical degrees of freedom, rather than being described by 
dynamical quantum fields (in particular, a quantized metric field). These fundamental ques¬ 
tions deserve further study with a view towards applications in cosmology. 


2 Expression for the electric current density in the presence 
of an inhomogeneous axial chemical potential 


2.1 States describing local thermal equilibrium 


In order to account for the local nature of the asymmetry between left-handed and right- 
handed particles, we propose to introduce a space- and time-dependent axial chemical po¬ 
tential, p^{x), and generalize Eq. (5) to a local relation between the pseudo-scalar density ps 
and /is, valid in the regime where /is <C T. We assume that local thermal equilibrium (LTE) 
is reached at length scales small as compared to the scale of spatial variations of /i 5 (x), i.e., 
^LTE ^ /W 5 (x)/|V/i 5 (x)|), and that the scale of variation of electromagnetic fields and matter 
flows is much larger than /lte> too- 

A quantum-mechanical system in local thermal equilibrium can be described similarly as 
one in perfect equilibrium; (see, e.g., [54]). Its state is tentatively described by the density 
matrix 


Plte = 2 ^ exp 




(9) 


Here l-L is the Hamiltonian of the system, and Z is the partition function (chosen such that 
tr(pLTE) = 1); is as in Eqs. (2) and (4), and T is the temperature (here assumed to be 
constant throughout space-time). 

The integration in the exponent on the right side of (9) extends over a spatial hyperplane at 
fixed time. Expression (9) may be generalized as follows: Eet (u^{x)) be the hydrodynamical 
4-velocity field of the plasma; (the condition that u^{x) = 1, m*(x) = 0, for i = 1, 2, 3 means 
that we are working in a co-moving coordinate frame). Eet (Ty,fix)) denote the components of 
the energy-momentum tensor of the system, and let /3(x) = {^kBT{x))~^ be a (possibly slowly 
space- and time-dependent) inverse-temperature field. Eet S be some space-like hypersurface, 
and let da^{x) denote the (dual of the) surface element at a point x G S. The covariant form 
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of the expression (9) is given by (see e.g. the book [54, Chap. 24]): 


Plte = ^ exp 


/ 

u^{x)Tp^{x) - y,5{x)J5u{x) 

/3{x)da’^{x) / 

s 


J 


( 10 ) 


An equivalent expression is given in Seet. 2.2, Eq. (32). Expeetation values of operators in the 
state introdueed in (10) are denoted by ((■))lte- We recall from Eq. (2) that 

(y 

which, by Eq. (1) (chiral anomaly), is conserved but not gauge-invariant. Under a gauge 
transformation, A,^ ^ Ai, — dyt], where r] is an arbitrary (smooth) function on space-time, the 
term 

J /3{x)fi5{x)J5^{x)da''{x) = J (3{x)ti^{x)J^{x)day{x) 

s s 

changes by 

j e''^P^d^{l3{x)^i^{x))Fpx{x)rj{x)da:,{x), 

as follows by integration by parts; using that d[pFpx^ = 0. Requiring that expression (10) 
be gauge-invariant (i.e., independent of the gauge function rj, for arbitrary rj), we find the 
constraint 

{/3{x)h5{x)) Fp^[daf^{x)]j^ = 0 , ( 11 ) 

or, in the language of differential forms. 


d(/3/i5) AFls^O. (12) 

If j3{x) = (3 is constant and S is a hyperplane at fixed time then (11) implies that 

B{x) ■ V/i 5 (a:) = 0, (13) 


i.e., gauge invariance implies that the gradient of the axial chemical potential is orthogonal 
to the magnetic induction. Condition (13) has a simple physical interpretation discussed in 
Sect. 2.3, below. 

Eet J^^{x) denote the electric (vector) current density, and let be as in Sect. 1.1; (see 
Eq. (1)). We set 

j^{x) = {J^{x))uTE, 35 i^) ~ ('^^(^))lte, (14) 

for yU = 0,1, 2, 3. In [1,2], the equal-time commutation relations between and Jg, 

(y. —* —* 

= -B{yA) -y), (15) 

71 

see [51], have been used to derive expression (7) for the B- dependence of the electric current 
density j: 

j{x) = -fX 5 B{x), (16) 

71 

assuming that does not depend on x and that aE = 0. We sketch how (16) is derived from 
(15) and then propose a generalization of (16) where /is may depend on space and time. 
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The conservation of the vector current density J^(x) implies that there exists a ’’current 
vector potential”, $(x), with 

J0(x) = -V ■ $(a;), J{x) = dt${x). (17) 

These expressions are invariant under the transformations $—)■$ + V A X, for an arbitrary 
vector field X independent of time. The commutation relations of the field with the axial 
charge density Jg are then given by 

= -B{y,t)5{x - y) + A n(f, f), (18) 

TT 

where the second term on the right side of (18) is an ’’integration constant”. Let vf denote 
the expectation value of n in the state ((■))lte of the plasma. If this state is homogeneous (in 
particular, (i?)LTE = 0) then 7f(x, y, t) = Tf{x — y, t) must be space-translation invariant. It is 
plausible that vf is independent of 5, so that vf is space-translation invariant even if (i?)LTE 7^ 0. 
In the following, we will write B for (i?)LTE, as above, in order to shorten our notation. Taking 
the expectation value of both sides of (18) in the state 

((■)) = ((■))lte 

of the plasma and integrating over y, we find that 

{[Q,Mx)]) = -B{x), (19) 

TT 

because the expectation of the second term on the right side of (18) integrates to 0 if vf = (If) 
is translation-invariant. 

Since $ commutes with the electric charge operator Q, eq. (19) implies that 

{[Ql/r,Hx)]) = ±^B{x), 

where Ql/r = |(Q ± Qs). 

Next, we observe that 

j(x) = ($(a:)) = i tr (plte[^, $(x)]) , (20 ) 

where Plte is given by (9), with y, 5 {x) = /is independent of x. 

Formally, the right side of (20) appears to vanish for a constant /is, because Jg is a con¬ 
served current. However, the field $ is so singular in the infrared that the formal calculation 
is deceptive. The right side of (20) must be regularized in the infrared by adding a mass term 
to the Hamiltonian, V. —> T-imass- Then Jg is not conserved, anymore, and \Hmass, Qb] 7^ 0. 
However Plte continues to commute with T-Lmass - t^b j Thus 


]{x) =(<F(x)) = itr Plte 


Ur 


t^b / fb{yU)d y + t^b / fb{yU)d y,^{x,t) 


a 


( 21 ) 


=/is([Qs,*^>(a;)]) = -y-bBix), 

71 


by Eq. (19). 
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Next, we eonsider the general situation where /is may depend on time and spaee. Let 
7 r{x, y,t) := (n(a;, y,t)), where H is as in Eq. (18). Then Eq. (21) generalizes to 

f{x) = ^Ij, 5 {x)B{x)+ A j Ti{x,y,t)y5{y,t)d^y, (22) 

Note that the eurrent j in (22) is conserved as a consequence of Eq. (13); (gauge-invariance of 
Plte)- If 1^5 is independent of spatial coordinates and the distribution vr is translation-invariant, 
the second term in Eq. (22) vanishes, and this equation reduces to (21). 


2.2 Axial chemical potential and axion field 

The second term in the expression (22) is non-local and the distribution 7 r(x, y, t) has not been 
defined, yet. We attempt to find a local expression for the four-current density j^. This current 
density should be proportional to (E)(i/) and linear in (/is); and it should be conserved, (i.e., 
satisfy the continuity equation). By introducting a pseudo-scalar field, 6 ^{x), indeed an axion 
field, related to the axial chemical potential y^{x), we arrive at the following ansatz for 

(23) 

with 05 proportional to /is in co-moving coordinates. The current defined by Eq. (23) 
is automatically conserved; (j'axion I® dual to the 3-form ^d^s A F, which is closed, because 
= 0 = dF). Eor to transform as a vector under parity and time reversal, 05 must be a 
pseudo-scalar field. 

To define its relation to /is let us separate the 0-components from the spatial components: 

Jaxion = - V05 ' B, jlxion = " U^B + V05 A e) . (24) 

Comparing these equations with Eq. (16) (see also Eq. (7)), we argue that ^^{x) = 05 (x). The 
identification of the axial chemical potential with an axion field ( 05 ) has first been proposed 
in [1,2]; (later it also appeared in the context of the quark-gluon plasma in Refs. [55-60], with 
the pseudo-scalar field being identified with a variable QCD theta angle). 

We wish to explain the relation between an inhomogeneous yfix) and 9^{x). Eq. (24) 
holds in any coordinate system, if we define Ei = Fq* and Bi = ^CijkF^^. However, in plasma 
physics, it is convenient to express F^^, in terms of electric and magnetic fields in co-moving 
coordinates: 

(25) 

where are 4-vectors of electric and magnetic fields orthogonal to the 4-velocity field 

u^, (i.e., = 0). Eor details see e.g. the review [61]). Taking F^^, in the form (25) 

and plugging it into the Eq. (23), we find that 


] axion 


a 

71 


{ 65 + v- 'V 6 fi)B{x) + (V 05 + v 6 fi) A S(x) 


+ 0{v^ 


(26) 


(where S,B axe the spatial components of the 4-vectors E^,B^. In what follows, we denote 
these fields again by E, B). Comparing Eq. (26) with (22) gives us the desired relation be¬ 
tween /i 5 and 05 

/i5 (f, t) = u'^df,95 (f, t ). (27) 
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We note that if the relation /i 5 (f) = 0^{t) holds in eo-moving eoordinates then Eq. (27) follows 
by transformation to the laboratory frame. 

Eqs. (27) and (26), (24) provide a loeal expression for the eleetrie four-eurrent density 
involving the axion field 6 * 5 . Our expression for explains why the electric current density 
i in (22) is non-local in the axial chemical potential /is. The current density j is odd under 
time reversal, T. The combination is T-odd, too, because B is T-odd and is T- 
even. The electric field (E) is T -even, hence a term proportional to (spatial derivatives of) /is 
and E would have the wrong transformation properties under time reversal. However, when 
introducing the axion field 6*s, it is easy to write a local combination of (spatial derivatives of) 
^s and E that transforms correctly under time reversal, namely E, which is odd under T, 
because the axion field is T -odd, (as the expression in (27) shows). Clearly, V^s is a non-local 
functional of /is. (Thus, introducing an axion field enables one to find a local expression for 
the electric current, instead of the non-local one on the right hand side of Eq. (22), and the 
second equation in (24) is the correct generalization of Eq. (21) to a spatially inhomogeneous 
plasma. 

An alternative route towards understanding the origin of the axion field 6 *s is the following 
one, (see [1, 2]): We consider a system of massless fermions coupled to the electromagnetic 
field in a state corresponding to a non-vanishing axial chemical potential /is. The action of 
this system is given by 


S[A,'il)\ 





(28) 


One can remove the term proportional to /is by performing a “chiral redefinition” of the 
fermion fields 

f (29) 

where 0^ = (see (27)). This transformation yields a Eujikawa Jacobian [62] 

J[A^, 0s] = exp (i (30) 

in the functional integral, yielding the modified action 

S|. 4 ,x,« 5 ] = jd‘x + + . (31) 


We see that if the axial chemical potential is spatially homogeneous then it can be re-expressed 
in terms of a pseudo-scalar field, 9^, that has an axion-like interaction with the electromagnetic 
field. 

Eq. (27) enables us to present the general expression for the density matrix describing 
local thermal equilibrium of an inhomogeneous plasma in an arbitrary coordinate system. 


Plte = Zj, ^ exp 


u^{x){T^fx) - d^9fx)J^fx)) /3{x)d 



(32) 


where is the 4-velocity field of the plasma and the expression plays the role of a lo¬ 

cal axial chemical potential jJi^{x). Expression (32) generalizes (10) to an arbitrary coordinate 
system. 

In Section 3 we will search for the correct dynamics of the axion field. But, before ad¬ 
dressing this problem, we propose to clarify the significance of the constraint (13). 
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2.3 Constraints on 6*5 

In this subsection, we tentatively interpret as a classical field, more precisely as a gener¬ 
alized thermodynamic parameter that labels states describing local thermal equilibrium with 
left-right (chiral) asymmetry. As discussed in Sects. 2.1 and 2.2, the requirement that the 
density matrix Plte, supposed to describe local thermal equilibrium (see (10), (32)) be gauge- 
invariant yields the constraint 


d{/3fi5) AF\^ = 0, (33) 

with 

(34) 

where u^{x) is the 4-velocity field of the plasma appearing in expressions (10) and (32) for 
Plte- Eq. (33) is identical to (11)-(12), and Eq. (34) follows from comparing expressions (10) 
(Sect. 2.1) and (32) (Sect. 2.2). In co-moving coordinates, vP{x) = 1, u{x) = 0. Choosing S 
to be a hyperplane at fixed time (in co-moving coordinates), and assuming that the temperature 
T is independent of x = (x, t), (33) becomes 

B{x) ■ V/i 5 (x) = B{x) ■ V9^{x) = 0, (35) 


which is Eq. (13) 

It may be somewhat surprising that the constraint (35) must be imposed on the choice of 
/is. In order to illustrate the origin of this constraint, we consider a plasma in a homogeneous 
magnetic field B ^ D. We choose an initial state where all left-handed (spin parallel to 
momentum) charged particles are located in the half-space x^ < 0 , while the right-handed 
ones (spin anti-parallel to momentum) are located at x^ > 0. We assume that the particles are 
non-interacting. Then the one-particle states are labelled by the quantum number of a Landau 
level and the component of the momentum parallel to B. If V/is is perpendicular to the field 
B then the densities of left-handed and right-handed particles are independent of time. But 
if .B ■ V/is 7 ^ 0 then left-handed and right-handed particles start to mix, as they move in the 
direction of B. As a result, the chiral charge density changes in time, which implies that local 
thermal equilibrium is lost; see Eigure 1 . 

As a second example, we choose /i 5 (x) to be non-zero and constant in the slab defined 
byO ^ x^ ^ L and /i 5 (x) = 0, elsewhere. We consider a B-field parallel to the x^axis. An 
electric current then flows from x^ = 0 to x^ = L, thus creating charge densities proportional 
to —B5{x^) and B5{x^ — L), respectively, on the two boundary planes of the slab. These 
surface charges create a non-vanishing electric field E{x) inside the slab, which is parallel to 
the x^-axis. Hence E{x) ■ B{x) 7 ^ 0, in the region where /r 5 (x) 7 ^ 0, i.e., inside the slab. By 
Eq. ( 6 ) (Sect. 1.1), l/xsl then decreases in time, with /r 5 (x, f) —)■ 0, as f —)■ cxd. Obviously, local 
thermal equilibrium is then only approached in the limit where t ^ 00 . 

2.4 Chiral asymmetry and electro-neutrality 

To conclude this section, we discuss consequences of electro-neutrality in a high-temperature 
plasma and in the presence of an axion field. Electro-neutrality means that the total electric 
charge density, jtot(x), vanishes. According to Eq. (24), the light charged particles (electrons 
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and positrons, with T 3> rrie) make a contribution 

Jaxion = ' B (36) 

71 

to the electric charge density. In a plasma containing only electrons and positrons, electro¬ 
neutrality implies the absence of matter-antimatter asymmetry, and Jaxion must vanish identi¬ 
cally. If = 0 then the 3-divergence of jaxion must vanish, i.e., 

0 = V ■ Jaxion = -(V4 ■ 5 + V 05 ■ 5) = -dt{Ve, ■ B), (37) 

71 71 

by (36). Expression (22) for the current density j implies that V ■ j = 0, as a consequence of 
Eq. (35). As a consequence, electro-neutrality, i.e., jaxion = 0 follows. 

If the plasma contains several types of charged particles (for example, and p, p) then 
electro-neutrality does not imply any symmetry between matter and antimatter, and there is 
no reason to expect that jaxion^ us given in (36), vanishes identically! One must then introduce 
electric 4-vector current densities for leptons and baryons and axial chemical potentials for 
all species of very light, charged particles. Eor simplicity, we consider a plasma consisting of 
electrons, positrons, protons and anti-protons only. Eet denote the electric current density 
of protons and anti-protons, and let 

i? = (38) 

Electro-neutrality then implies that 

= -jJxion = — V05 ■ i? (39) 

TT 

Eet be the total lepton vector current density, and let denote the total baryon vector 
current density. These vector currents are separately conserved; hence 

J^B - Jl (40) 

is conserved. In order to describe the matter-anti-matter asymmetry observed in the universe, 
we introduce a chemical potential /Tb-l = ^b-l conjugate io Ql — Qb, where 

Qlib= j Jl/Bix,t)d^x (41) 

t=const 

are the conserved lepton and baryon electric charges. If the distribution of asymmetry between 
matter and anti-matter in the universe turned out to be inhomogeneous the scalar quantity 6 *b-l 
would have to be taken to depend on time and space, and a term 

J {Jb{x) - Jl{x)) (3{x)dau{x) (42) 

s 

must be added in the exponent of the expression for the density matrix that describes local 
thermal equilibrium, on the right side of Eq. (32); (see Sect. 2.2). To complete the picture 
we will have to find an equation of motion for the scalar field ^b-l- We note that there are no 
constraints on 6 'b-l similar to (33), (35); but there is an equation analogous to (5) discussed in 
the next section. Before we turn our attention to the dynamics of 6 'b-l, we study the dynamics 
of 05 . 
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Figure 1: The simplest eonfiguration with a non-zero gradient of The region < 0 

has only left-handed partieles and the region > 0 has only right-handed partieles. If we 
consider a configuration where the magnetic field B has a non-vanishing 1-component, Bi, 
the left-handed particles drift to the right, while the right-handed particles drift to the left, 
(because their state is parametrized by a non-zero pi). Thus, the state is not on equilibrium 
state. But if B is perpendicular to the x-axis the particles move parallel to the boundary plane 
{x = 0}, and the state is an equilibrium state. 


3 Search for the dynamics of the axion field 

Having introduced a time-dependent axial chemical potential it is necessary to specify its 
dynamics. Eqs. (23) and (24) suggest that it may be convenient to derive the dynamics of /is 
from the one of the axion field 9^. 

In the homogeneous case, i.e., when = 65 ( 1 ), Eq. (1) can be rewritten as 

4:71 

where A is a constant with the dimension of an ‘energy’; (the field 6*5 is dimensionless, because 
/is = 0s must be an ‘energy’). By Eqs. (1-6), we have that 

T‘2 

A' = TT’ (44) 

where T is the temperature of the plasma. (We have used that = 8.E ■ B, and 

(■) indicates spatial averaging.) In (43), {F^^) is the expectation value of the electromagnetic 
field tensor in the state ((■))lte- We should ask what the correct equation of motion of 0s is 
that, upon taking spatial averages, yields Eq. (43). The axion field 0s has been introduced in 
order to describe plasmas exhibiting a chiral asymmetry; in other words, 0s may label states 
of a plasma - as do other thermodynamic parameters, in particular the temperature. Eor this 
reason, the equation of motion of 0s has no reason to be relativistic; (i.e., to preserve its 
form under arbitrary coordinate transformations). It may therefore be plausible to consider a 
diffusion equation for 05 , which we discuss below. 

The chiral anomaly tells us that 

TT 
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see Eq. (1), Seet 1.1. Taking the expeetation in the state of the plasma (and viewing E and B 
as elassieal fields) yields 

=—E-B, (45) 

TT 

with = {J 5 )T,ti 5 - On the basis of Eqs. (43) and (45) one argues that 

df,j5 = (^with , (46) 

where e is a term whose spatial average vanishes, i.e., e[6^] = 0 . 

If 6*5 is interpreted as a kind of thermodynamie parameter labeling an inhomogeneous state 
of the plasma in loeal thermal equilibrium then the following ansatz for jg as a funetion of 6*5 
in eo-moving eoordinates is reasonable: 



3l = k%, Jg = A2DV(0g), 

(47) 

and 






A2 = 

'931 

, with 75^5=0 = 0 . 

M5=0 

(48) 

Note that these equations imply that 






DV3l=35, 

(49) 


a relation that ean also be eonfirmed by direet eomputations; (see e.g. [63]). In eo-moving 
coordinates, (45) and (47) then yield the equations of motion 

A2 (e, - 

which imply (46). Together with Eq. (27) we find that 

2a ^ ^ 

A^ (/is - DA/is) = —^ ■ S, (51) 

TT 

in accordance with Eq. (6), which is an inhomogeneous diffusion equation for jj,^. This diffu¬ 
sion equation can also be obtained from a more foxmdX first order expansion of the constitutive 
equation for the global current jg ; (see e.g. [64, 65], or the book [66]). Taking into account 
chirality flips due to small masses of the charged particles, Eq. (51) should be generalized to 
the equation 

(As ~ -DA/ig) -f Ej/ig = — E ■ B. (52) 

71 

Eurthermore, local thermal equilibrium imposes the constraint 

V/i5 ■ 5 = 0, (53) 

see (13). The electric vector current density is given by 

Jtot = Jaxion + 7p = - V05 ' B + = 0, (54) 

^71 ^ 

jtot Jaxion “1“ jp Jaxion “1“ johm “ ^ (jE^ (55) 
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with \jp\ \ji\, because protons are much heavier than electrons and positrons. 

Of course, we also require the Maxwell equations, with 3^ = jL given by (54) and (55). 

Assuming that the baryonic current jp can be incorporated in johm = o'E, as we have done 
in (55), we have arrived at a complete system of field equations for 6 * 5 , E and B\ 

The conductivity (tensor) a must be calculated from transport equations for the plasma. 
Apparently, the equation of motion for the scalar field ^b-l that determines the matter-anti- 
matter asymmetry is not needed here; although it is of considerable interest interest to derive 
one, (as will be discussed elsewhere). 

Equations (52)-(55) are supposed to be valid in co-moving coordinates. In general coordi¬ 
nate systems equations (53) and (54) must be replaced by 


d5 = K.0^ + u- V65 

(56) 

and 

JL + U-Jtot = 0; 

(57) 

see Eq. (27). 

The analogue of Eq (5) for the field 6 *b.l is the relation 


3% - 3l = 

(58) 


where (■) denotes spatial averaging, and e: is a constant with the dimension of an energy. Since 
3b ~ 3l ^ conserved current, it follows that 


0B-L = 0 (59) 

A relativistic equation compatible with (59) is 

□0B-L = 0, (60) 

which describes a free, massless scalar field. As in our discussion of the equation of motion 
of 6 * 5 , one could also envisage an equation of the form 

Ab-l ~ f9A/XB-L = 0 , (61) 

with /ig-L = ^B-L in co-moving coordinates. Choosing the hyperplane S appearing in Eq. (42) 
to correspond to a constant time t (in co-moving coordinates, with = l,u = 0 ) and impos¬ 
ing the condition that expression (42) is independent of t (i.e., the matter-anti-matter asymme¬ 
try does not change in time, anymore, in the old universe), we find the constraint 

(- Jf) = 0. (62) 

3.1 Dynamics of a fundamental axion field 

It is advisable to also consider the possibility that 6*5 is a fundamental field evolving in time 
according to some relativistic Hamiltonian dynamics. The equations of motion for 6*5 can then 
be derived from an action principle, the action of 6*5 being denoted by 

Seff[A,e,], (63) 
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where A is the electromagnetic 4-vector potential. In order for the action (63) to reproduce 
the desired expressions for given in (23-24) (Sect. 2.2) , it must have the form 


a 


S,ff[A,e5] = ^o[05] + ^J 

The equation of motion for 6 ^ derived from (64) is given by 


per- 


69 , 


= 0 , 


or 


69 , 

By (45), this implies that 


6Seff[A,95] 6 So[ 95] ^ p 


69, 271 

6S,ff[A,9,] 5.So[05] 


69, 

A reasonable ansatz for Seff[A, 9,] is 


69, 


= 


71 


Seff[A,9,] = I d^x [A\dp9,){x){d^9,){x) + U{9,{x)) - ^9,{x)e^’'^PFp,{x)F^p{x) 


(64) 


(65) 

( 66 ) 

(67) 


This action yields the equations of motion 


cy 

A^n9, + U'{9,) = —e^^'^P-Fp^Fp^ . 


( 68 ) 

(69) 


If 9, is spatially homogeneous and U = 0 this field equation reduces to (43). Eq. (69) has the 
form (46), provided U'{9,) = 0; (e.g., if (/ = 0). Eq. (69), with a term T f9, added on the left 
side, represents an alternative to Eq. (52). 

We conclude that, in a system with a fundamental axion field, the chiral MHD description 
of the plasma can be expected to be applicable even if the fermions are not ultra-relativistic. 
Eor example, if an axion-like particle exists and represents a contribution to dark matter, then 
the production of dark matter in such a system may be accompanied by the generation of 
magnetic fields (see for example [67]). The phenomenological viability of this scenario will 
be discussed elsewhere. 

If Eq. (69) is the correct equation of motion it is temping to argue that not only F^y cor¬ 
responds to a quantized field, but the axion field 9, must be quantized, too! Eurthermore, the 
action (68), with U = 0, suggests to unify the electromagnetic gauge potential A and 9, to a 
gauge potential, A = {A, A9,), on a slab of width oc = const m five-dimensional space- 
time, with the inverse temperature; (see (44)). Then the term f d^x {A^{dp9,){x){dP9,){x)} 
originates from a contribution to the five-dimensional Maxwell term, after dimensional reduc¬ 
tion. These ideas are discussed in more detail in Appendix A . 


4 Dynamics of chiral MHD in the inhomogeneous case 

It was demonstrated in [3] that the presence of a dynamical field in the Maxwell-axion 
equations leads to the generation of an “inverse cascade” - i.e., to the transport of magnetic 
energy and helicity from shorter to longer scales. In what follows we present a preliminary 
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analysis of the inhomogeneous equations derived in Seet. 3 and show that the behaviour of 
solutions is qualitatively similar to that of spatially homogeneous solutions, as diseussed in [3]. 
In our analysis, we negleet 

1. the diffusion term A/xs in Eq. (52), assuming that the timeseales are much longer than 
typical diffusion times; 

2. the coupling to the velocity field u of the plasma, as in Eqs. (56)-(57); i.e., we set m 0; 

we also neglect, as usual. Maxwell’s displacement current proportional to E; 

3. the chirality flipping term in Eq. (52), i.e., we set E/ = 0. 

We thus analyze the following simplified system of equations: 

VAB = crB + ^(flsB + vflj X i), 

A% = —E-B, ( 70 ) 

TT 

5 ■ V05 = 0, 


with = 65 . 

4.1 Exact single-mode solution 

In this section, we set out to find the simplest (but non-trivial) exact solutions of Eqs. (70). 
These solutions can be considered to be direct generalisations of the “single-mode” solutions, 
considered in Ref. [3], to the inhomogeneous case. 

We start with the following ansatz for the magnetic induction B: 

B{x) = 5(f) (sin(fcz), cos(/c 2 :), 0 ) (71) 

Eor this choice of 5, we have that 

W AB = kB, (72) 

which is a special case of the-force-free configuration first described in [ 68 ]; (see also [69,70]). 
A general force-free configuration (72) is given by the following expression [70]: 

B = k~^'V (efi) -f V A (efi) (73) 


where e is an arbitrary unit vector and is a. solution of the Helmholtz equation: 

Afi + fcV = 0. (74) 

The configuration (71) corresponds to the choice e = (1, 0,0) and = sin(/c 2 :). Therefore, a 
tracking solution, found in [3], 

{ ^5 = -k 

E=0 ^ > (75) 

B =5o(sin(/c2;), cos(fcz), 0 ) 

is also a solution of the inhomogeneous Maxwell equations (10) for any value of the constant 
Bq. However, if one starts with non-zero 5(0) or with /X 5 7 ^ ^k (or if the initial conditions are 
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not monochromatic, see the next subsection) the value B^, to which the solution B{t) arrives 
at late times, will depends non-trivially on the initial conditions as we demonstrate below. 

By the Bianchi identity, the electric field corresponding to the configuration (71) is given 
by 

E = cos{kz), 0), (76) 

k 


hence 


E{x) ■ B{x) 


2 k dt 


B\t) 


is space-independent. 


(77) 


Next, we plug this ansatz into the system of Eqs. (70). This yields 


« dB‘^{t) 
dt 2Tik dt 


(78) 


or 


/^5(t) = /i5 


a B^{t) - Bl 
271 kA? 


(79) 


where Bq = i?(0) and /Xg = /i 5 ( 0 ). This is the exact solution for ^^{t), for all times t, derived 
from the one-mode ansatz (71). Such a solution has been first described in [3], and we see 
here that the local nature of the axial anomaly equation (70) does not invalidate the conclusion 
in [3]. 

Eqs. (70) allow us to determine B{t) and its asymptotic behaviour for large t. This is seen 
by plugging expression (79) for /i 5 (f) into the Maxwell equations, which then yield a closed 
system of equations for B: 


V AB 


A “ 
aE H— 

71 



c^ B\t)-Bi \^ 

71 2kA‘^ J 


(80) 


Using equation (76), and recalling that V AB = kB, for B as in (71), we arrive at a non-linear 
ODE for one scalar function B{t)\ 


kB{t) 


O' A/ X Oi 




a B‘^{t) 


Bl 


7T 


2 kA^ 


B{t) 


( 81 ) 


This equation can be integrated explicitly and yields B{f) as a function of Bq and /ig and we 
can write down the explicit time-dependent solution of the system (70): 


Bit) 

< 


Cl 


^1 + Uaexp 



a B\t) - Bl 

271 


where we have introduced a parameter 7 given by 



(y j Q 

—kfx^ + 
71 



with 


0-2 R2 
2 _ -Dq 

PO — ^2A2 • 


(82) 


(83) 


The constants Ui,U 2 can be expressed in terms of i?o,/^ 5 ,and k, we will not provide their 
explicit form. 
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Next, we analyse asymptotie behaviour of B{t) as t ^ oo. To this end we put B{t) = 0 
and find an algebraie equation for B^o = B{t —)■ cx)): 


^Bqq 


a 

TT 



a Bl,-Bl \ 

TT Ikh? ) 


B„ 


(84) 


There are two distinet solutions of Eq. (84). The trivial one, i.e., i?oo = 0, is approached iff 


Trivial solution ^ > 72 


or 


k > 


s • \/f * (#)■ 


(85) 


Notice, that although B^a = 0, the —)■ 00 ) can still be non-trivial. In the case of zero 
magnetic field the asymptotic form of the /is (f) is not bound to be equal to ^ and its value is 
determined by the Eq. (79). 

The non-trivial solution of Eq. ( 86 ) exists if fc < 7 . It is given by 



Bi 


1 + 




where /i“ = —>■ 00 ) 


Trk 

a 


( 86 ) 


The necessary condition for the non-trivial solution ( 86 ) to exist, i.e. > 0, is given by 



(87) 


It is exactly opposite to the condition (85), as it should be. 

In summary, we were able to find an exact solution of the non-linear system of differential 
equations (70) and found condition (85) (or equivalently (87)) that relates the initial value of 
= 0), B(t = 0) and the wave-number k and determines the ultimate fate of the single¬ 
mode solution. 

Notice that if, initially, = ^k = /i“ then the solution is stationary, i.e., Bq = Boo- 
If, however, ^/Xg 7 ^ k then b{t) will increase or decrease and hence cause ^^{t) to change in 
time, as is seen from Eq. (79) . Indeed, if 7/^5 > k then h{t) > 0, and the amplitude of the 
.B-field grows, while, according to (79), ^^{t) decreases - as it should. But if 7/^5 < k then 
b{t) < 0, i.e., the amplitude of the 5-field decreases, and, according to (79), /i 5 (f) increases 
towards /i“. To demonstrate this, we consider the behaviour of the solution for small times 
when B(t) ~ Bq and ~ /Xg. To this end we linearize the system of Eq. (81) around some 
special solutions. We set B{t) = Bq{ 1 + eb{t)) where 6(0) = 0. To first order in e we obtain 
the following equation for 6 (f): 


a 


-dl-k] = 


TT 


- k) 6(f) - ^6(f) - ^6'(f) 


+ 0{e). 


The solution of Eq. ( 88 ) is given by 


1 — exp 




6(f) = 






. A OilA 

where A/r = —^ — k. 


2kAfi 


TT 


and it demonstrates the conclusion. 


( 88 ) 


( 89 ) 
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4.2 Stability analysis of the helical single-mode solution 


In this section we sketch an analysis of stability of the solution found in the previous section. 
To this end, we assume that the initial configuration has reached its time-independent form 
with some k, B^o and /i“. We then perturb the vector potential A and the chiral chemical 
potential /xs by 


/i5(T,t) 

A(x, t) 


=-h dn5{x,t) 

a 


Bo{z) 


+ a{x, t) 


(90) 


where Bo{z) has the form (71), with amplitude B{t) = Boo given by Eq. (86). In (90), 
fi 5 {x, t) = /r“ if 5 fi 5 {x, f) = 0; see (86). 

The (equation expressing the) chiral anomaly (see Eq. (70)), linearized around the back¬ 
ground /r 5 (a;) = a~^TTk and A{x) = k~^Bo{z), then yields 


which can be integrated to 


dt 


2 a 

TT 


Boiz) 


da 


6fi5{x,t) 


2 a 1 

71 


Bo{z) ■ (a(f, t) - do{x)) + 5/i5(^, 0) 


(91) 


(92) 


where do{x) := a(x, f=0). Below, we choose as initial condition 6 fi 5 {x,t = 0) = 0. This 
is obviously a special choice. However, physically, it is the most interesting one, because the 
spatial inhomogeneity of the chiral chemical potential is then caused by fluctuations of the 
magnetic field. The analysis of a more general situation will follow. Using (92), the equation 
of motion for a is seen to be given by 

- Ad= + kV A (a) - "^^(^Boiz) ■ d - Bo{z) ■ do^Bo{z). (93) 

We analyse the special solution found by assuming that a depends only on the spatial coordi¬ 
nate X and on time t; i.e., a(x, t) = [ax{z, t), ay{z, f), O). The equations then reduce to 


a'^ = adx + ka'y + 2l3‘^{AT?{kz){ax — a°) -f sin(fcz) cos{kz){ay — a^)), 
a'y = ady — ka'^ -f 2(3“^{sm{kz) cos{kz){ax — a°) -f {kz){ay — a°)), 

where (3'^ is given by 

.2 _ o^^Bl 

^ ~ 7r2A2 • 

(similar to the definition of /Sq in (83), but with Bq replaced with Boo)- 
We rewrite the system of equations (94) as a matrix equation 


(94) 


(95) 


fadt + l3‘^-dl kd^ \ fax\ ^ _ r .2 (- (^os{ 2 kz) sm.{ 2 kz)\ 

y —kdz adtAl3‘^ — d‘lj \ayj y sm{ 2 kz) cos( 2 kz)J yay — ayj’ 

' -V-^ -V-' 

=L = 0 (z) 

(96) 
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or, schematically, as La = —(3‘^0{z){d — do). Defining U by 


U = 




( 97 ) 


we find that 


U-^LU 


and 


Furthermore, 


/ adt + (3"^ — dl — ikdz 0 

y 0 adt + — ^2 + ikdz 

jy—1 {^ fO^x ^Ojy\ f 

\ay) 2 \a^ - iayj \a-J 


U-^ 0 {z)U 


0 

^ 2 ikz 


^—2ikz 

0 



(98) 


(99) 

( 100 ) 


We then have that 

L+a+ = /32e-2*^"(a_ -a°), 

. ( 101 ) 

L_a_ = ^(a+ — a° ), 

and we choose as an initial condition a plane wave: a° = a° = hence VAuq = Qdo- 
One may attempt to solve Eqs. (101) by iteration. In zeroth approximation, we set the right 
side of Eqs. (101) to zero and get 


L+af = 0; af (0,.2) = e-*'?", 

L_a® = 0; a^°\0, z) = 

These equations are solved by 

af = a™ = A = 

a 

The condition that perturbation grows in time, i.e. A > 0, is translated into 


A > 0 


k- Jk^- 4/?2 k+Jk'^- 4/?2 

--< g <-- 


( 102 ) 


(103) 


(104) 


That is, for (3'^ < fc^/4 and q obeying the inequalities in ( 104), a perturbation with wavenumber 
q grows in time; (i.e., the single-mode solution of Section 4 is unstable). We stress that, com¬ 
pared to the results of Ref. [3], not every mode, longer than 27i/k, can grow. The difference 
is given by the presence of term, which is a non-linear contribution from the background 
field. Boo- Eor (3^ > fc^/4, the solution of the previous Section 4.1 is stable with respect to the 
perturbation in the form 

To point out an important difference between the stability analysis of Eqs. (70) and the 
results of [3], we iterate Eqs. (101) again, yielding 


- 1); ai!^(0, z) = 0, 
- 1); aL^^(0, z) = 0, 


(105) 
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which is solved by (we only write the explicit expression for a_ being its complex con¬ 
jugate) 



p2^-i{2k-q)z 

2k{k — q) 





= (106) 


CJ 


for q < k. We see that is a wave with inverse wavelength 2k — q. If q < k (i.e., at 
f = 0 , the perturbation has a longer wavelength than the original single-mode solution) then 
2k — q > k\ i.e., a wave with shorter wavelength is excited. This result is a consequence of 
the non-linear nature of Eqs. (70) and the coupling between the modes (i.e. the generation of 
mode 2k — q starting from two modes with the wave-number k and q) is expected. 

Below we show that the coefficient in front of the term in is always smaller than 
the term in front of the harmonic. The ratio of coefficients can be found to be (we take 
times large enough so that \it 3> 1 and \t S> 1 ) 


a+\t,z)\ (3^ ^ q 

a+^(f, ;2)| Ait>i;At>i 2k{k-q) 2k 2 


(107) 


(in view of the condition that [3'^ < kq — that follows from A > 0, Eq. (104) and considering 
only modes with q < k). As a result at most | of the energy could be transformed into the 
short wavelength mode. 


Summary: Starting from the simplified system of equations of chiral electrodynamics, see 
Eqs. (70), we have shown that the helical single-mode solution is stable for sufficiently strong 
magnetic fields. Eor weak magnetic fields, the long-wavelength perturbation grows in time. 
A short-wavelength mode also gets excited, but its amplitude is always parametrically smaller 
than the one of the long-wavelength mode. Thus, transfer of energy and helicity from a sin¬ 
gle helical field mode to a long-wavelength perturbation, as described in [3], is observed for 
solutions of the full system of equations with an inhomogeneous axial chemical potential. 


5 Conclusion and outlook 

In a relativistic plasma of charged particles, the axial anomaly is accompanied by the appear¬ 
ance of new degrees of freedom described by a pseudo-scalar field 6 * 5 ( 0 ;, f). This field has 
the property that its time derivative is equal to the axial chemical potential /is parameterizing 
states of local equilibrium of the plasma. When coupled to the electromagnetic field and to 
the motion of the fluid the field 6*5 appears in a system of equations that we call the “equations 
of chiral magnetohydrodynamics”; (see Eqs. (47) through (54)). If 6*5 is spatially uniform and 
depends only on time 9^ is nothing but the axial chemical potential of the light charged 
fermions. In this paper, we have studied states of the plasma corresponding to local equilib¬ 
rium, with 6*5 depending not only on time t, but also on x. Our analysis is significant, because 
the presence of a non-zero (homogeneous) axial chemical potential /is leads to an instability of 
some solutions of the Maxwell equations accompanied by the generation of helical magnetic 
fields [1, 2, 5]. This instability leads to the appearance of strong electromagnetic fields that 
will, in turn, back-react on /is, making it spatially inhomogeneous. 

We have analyzed special solutions of the system of chiral MHD equations and have shown 
that the qualitative conclusions reached in [3] remain valid. In particular, we ave shown that 
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there is a stationary solution (helieal single-mode solution) of our system of non-linear equa¬ 
tions with the property that the growth eaused by a non-vanishing axial ehemieal potential 
exaetly eompensates the Ohmie dissipation (Seetion 4). Another important property of so¬ 
lutions in the presenee of a homogeneous axial ehemieal potential is the “inverse eascade” 
phenomenon, i.e., the transfer of energy and magnetic helicity from short to large scales, de¬ 
scribed in Subsection 4.1. This transfer is not caused by turbulence in the plasma, but is an 
effect derived from the chiral anomaly and caused by the chiral imbalance in the plasma . We 
have shown that, under suitable conditions (Section 4.2), this phenomenon appears to persist 
for spatially inhomogeneous fields. 

Our analysis can be expected to be relevant for the study of various physical systems 
exhibiting a chiral (left-right) imbalance, including the quark-gluon plasma and the plasma 
in the early Universe. In an analysis of the plasma in the early Universe, it might not be 
legitimate to neglect the coupling of the axion field 6*5 to the space-time curvature, (namely to 
the term Rij^uarRap^^ where R/M^ar is the Riemann tensor). We plan to return to this 

topic elsewhere. 

Many elements of our analysis remain unchanged when the axion is treated as a funda¬ 
mental field with relativistic hamiltonian dynamics. As discussed in Sect. 4, it can induce a 
left-right asymmetry in the distribution of electrically charged fermions and trigger the growth 
of magnetic fields. 
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Figure 2: Five-dimensional geometry. Light chiral fermions are localized on the boundaries 
of the five-dimensional slab {x'^ = 0 and x^ = t). 


A Electrodynamics on a slab in five-dimensional Minkowski 
space 


The connection between the chiral chemical potential and the axion field 6*5 becomes 
manifest if one studies electrodynamics on a slab M3.1 X X of five-dimensional Minkowski 
space, where X is a finite interval extending into the fifth dimension; (cf. [1, 2, 71]). Here 
we briefly review some important elements of this story. We consider very heavy charged 
four-component Dirac fermions populating the 5-dimensional “bulk”. After coupling these 
fermions to an external electromagnetic 5-vector potential A and, subsequently, integrating 
them out, we obtain a low-energy bulk effective action for the gauge field A and a bound¬ 
ary effective action for boundary degrees of freedom, which turn out to be massless chiral 
fermions localized on the (3 -f 1)-dimensional “top” and “bottom” boundary components of 
the slab (see Fig. 2)f The effective action is given by 


^(5) 


d^x 


— —-L ^ , 

4f 327r 


.abode A IP IP 
^aXbcXde 


boundary action for chiral fermions 
higher-order non-local terms, 


(108) 


where ^ is the length of the interval X, a is the 4-dimensional fine-structure constant, and Fab is 
the 5-dimensional field strength (a, 6 = 0,1, 2, 3,4). Fermions of opposite chirality are located 
on opposite boundary components, i.e., at = 0 and x'^ = i, respectively. We turn on an 
electric field pointing into the fifth dimension, i.e., F^a = 0, for a = 1, 2, 3, and X 04 = const.. 
Then there are non-zero charge densities, Pl{x) and pr{x), of chiral fermions located on 
opposite boundary components of the slab. These left- and right-chiral fermions couple to 
the electromagnetic gauge field restricted to the boundary hyperplanes. They give rise to 
anomalous currents on both boundary hyperplanes. The action in (108) must therefore 
contain the Chern-Simons term that cancels the gauge anomalies of the chiral fermions on the 
two boundary components. The bulk electric current corresponding to the action is given 
by 


■^These chiral fermions may acquire a mass through tunneling between the two boundary components, which 
shows that they may be used for purposes of reasonably realistic model building. We will ignore this possibility 
in the following. 
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It has a contribution corresponding to the Chern-Simons term 
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Note that is divergence-free in the five-dimensional bulk. Its divergenee does, however, 
not vanish on the boundaries, where it is proportional to respeetively, and 

eaneels the divergenee of the anomalous eurrents of the massless ehiral fermions loeated on 
the boundary eomponents of the slab - a phenomenon known as anomaly inflow. 

In order to reveal the eonneetion between this model of 5D eleetrodynamies and axion 
electrodynamies, we study the dimensional reduetion of the five-dimensional theory. We de¬ 
fine a sealar field depending only on the eoordinates X (yX t ij X ^ X ^ X of a spaee-time 

point X G X X by setting 


eflx) 


dx^ AflX) 
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where 7 is a straight path in the fifth direetion eonneeting the two points Xi = (x, x"^ = 0) and 
Xu = (x, x"^ = t) located on opposite boundary eomponents; see Fig. 2. We assume that the 
components A^, gi = 0,1,2, 3 of the gauge field A are independent of the -coordinate (or 
are averaged over x^). By Eq.(l 10) 


O5 = — dx'^ = / dx^E^. (Ill) 

J 'y J p 

The right side of this equation is the voltage drop between the two boundary eomponents, 
which is nothing but the difference of the ehemieal potentials of left- and right-ehiral fermions. 
The dimensional reduction of the action (108) yields exactly the effeetive aetion given in 
Eq. (68) -f a 4D Maxwell term. This implies that the axion-photon eoupling eonstant ean 
be identified with the size of the interval X extending into the fifth dimension. Comparing 
the equations derived in this appendix with Eqs. (5) and (6) of Sect. 1, we observe that the 
temperature T can be identified with the inverse of the width ^ of the slab. 

The theory outlined here ean be viewed as a five-dimensional eousin of the quantum Hall 
effeet. 


B Linear analysis 


B.l Growth of the long-wavelength modes 

Eet us now perturb the stationary solution for the one-mode (80) by a small long-wavelength 
mode Bp. 

Bflzfl) = eBobi{t)[sm{kiz),cos{kiz),0) (112) 

The linearised equation does not have the simple form Eq. (88) anymore. However the prop¬ 
erty (77) still holds for the mode Bi separately. What is not true, however is that the overall 
P 5 {t) is now spaee-independent. It is governed by the following equation (where again we 
introdueed B{t) = Bo(l + e6(f)) 
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To the first order in e Eq. (113) has the following form 
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We ean integrate Eq. (114) over time to get: 
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It is important to notiee that in this order in e 


V05OC (0,0,1) (116) 

— points in the ^-direction and therefore B -VO^ is indeed equal to zero. 

If /xs is homogeneous in spaee (ie. if we negleet the last term in Eq. (115)) we end us with 
the usual two-mode equation: 


kB{t) = -jB{t) + -fi5{t)B{t) 

k IT 

kibi{t) = + -fi5{t)bi{t) 

ki IT 


(117) 


If instead of solving equations (117) direetly we would expand them in e (in view of the 
subsequent non-homogeneous case) we would get the following. 


We start with the 


,,0 _ nk 

E5 - a 


(i.e. the one that would make a single mode stationary). We 


then see that the correction to this mode obeys an ODE (compare the last two terms in Eq. (88)) 


ab\t) =-^^b{t) (118) 

This equation has the solution b{t) = 0 if one starts from 6(0) = 0. This probably means that 
B{t) does not change in the linear order in e. In the first order in e the eq. for bi{t) has the 
form, expected from (117): 

<jb[{t) = {kki — k\)bi{t) (119) 

which gives exactly the exponential solution with 1^5 = ^- 

Qualitatively one has constant amplitude for the short wavelength solution (the non-trivial 
evolution appears at the order e^). The growth of the mode bi{t) starts immediately and goes 
on until the neglected terms become important. 


B.2 Two-mode inhomogeneous solution 

Eet us now turn to the case of two-mode inhomogeneous solution. The expression for /xs is 
given by Eq. (115) and the following equations for the different modes appear: 

The mode with the wave vector k {short wavelength) 


ab'{t) 


Bl 

^7! 


b{t) 


(identical to (118)) 


( 120 ) 
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The mode with the wave veetor ki {long wavelength): 


ah\{t) = {kki - kl)hi{t) + (hi{t) - 6i(0) j (121) 

Due to the non-linearity of the equations, we also have the very short wavelength mode {2k — 
ki) excited. If we assume that its initial amplitude was zero, we have 


0-62 (f) 


a‘^Bl{2k - ki) 

27r2fci/2 




( 122 ) 


— a term, similar to the last term in Eq. (121) but with an additional “enhancement” k/ki S> 1. 
However, this term is sourced by (&i(f) — &i(0)) and will not get excited until the mode hi had 
grown significantly. 
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